Dynamical Formation of the Dark Molecular Hydrogen Clouds 

around Diffuse H II Regions 

Takashi Hosokawa 1 and Shu-ichiro Inutsuka 2 



ABSTRACT 

q | We examine the triggering process of molecular cloud formation around dif- 

fuse H II regions. We calculate the time evolution of the shell as well as of the 
H II region in a two-phase neutral medium, solving the UV and FUV radiative 

<^ transfer, the thermal and chemical processes in the time-dependent hydrodynam- 

ics code. In the cold neutral medium, the ambient gas is swept up in the cold 
\ (T ~ 30 — 40 K) and dense (n ~ 10 3 cm' 3 ) shell around the HII region. In the 

shell, H 2 molecules are formed from the swept-up H I gas, but CO molecules are 

00 ' 

hardly formed. This is due to the different efficiencies of the self-shielding effects 
between H2 and CO molecules. The reformation of H2 molecules is more efficient 
with a higher-mass central star. The physical and chemical properties of gas in 



O 



the shell are just intermediate between those of the neutral medium and molecu- 
lar clouds observed by CO emissions. The dense shell with cold HI/H2 gas easily 
Q-f becomes gravitationally unstable, and breaks up into small clouds. The cooling 

layer just behind the shock front also suffers from thermal instability, and will 
fragment into cloudlets with some translational motions. We suggest that the 
predicted cold "dark" HI/H 2 gas should be detected as the H I self-absorption 
(HISA) feature. We have sought such features in recent observational data, and 
found shell-like HISA features around the giant H II regions, W4 and W5. The 
shell-like HISA feature shows good spatial correlation with dust emission, but 
poor correlation with CO emission. Our quantitative analysis shows that the 
HISA cloud can be as cold as T ~ a few x 10 K. In the warm neutral medium, 
on the other hand, the expanding diffuse H II region is much simpler owing to 
the small pressure excess. The UV photons only ionize the neutral medium and 
produce a warm ionized medium. 

Subject headings: Circumstellar matter H II regions - ISM: molecules - STARS 
: formation 
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Introduction 



The UV [hv > 13.6 eV) and FUV (11 eV < hv < 13.6 eV) radiations from massive 
stars are important for the cycle of the interstellar medium (ISM). Stellar UV radiation ion- 
izes the neutral medium, and this photoionization is the main heating process in the ionized 
medium . The diffuse warm ionized medium (WIM) is widely distributed throughout our 
galaxy (IHaffner et al.l 120031). and the phot o ionization is a prom ising process to supply the 
ubiquitous WIM (e.g., iMiller & Coxl Il993l : bove & Shul3 Il994h . The FUV radiation pho- 
todissociates the molecular medium, and photoelectric heating via the dust absorption of 
FUV photons is the primary heating process in the neutral medium. Since the chemical 
and thermal balances determine the equilibrium of the neutral medium for a given density, 
the strength of th e FUV radiation field iney itably determines the physical state of the neu- 
tral medium (e.g.. IWolfire et al.lll995l . 120031 ). The well-known equilibrium n-P curve shows 



that the two different phases can be achieved at some pressure levels. It has been gener- 
ally accepted that the observed cold neutral medium (CNM) and warm neutral medium 
(WNM) correspond to these two thermally stable states (see I Coxl 120051 . for a recent re- 
view). Besides its significance in the diffuse ISM, the UV radiation also works to form the 
cold and dense molecular structure. The shock front (SF) sometimes emerges ahead of the 
ionization front (IF) owing to high pressure in the H II region. The compression of the 
neutr al medium behind the SF is one of the possible pro cesses to form the molecular cloud 



Koyama fc Inutsukall2000l . 120021 ; iBergin et al.ll2004l ). The ambient medium is swept up 



(e.g. 

by the expanding SF, and molecules reform in the dense shell. The density a nd strength 



of FU V radiation in the shell determine the formation efficiency of molecules (jElmegreen 



19931 ). Some observat ions have shown clear features of this mode of triggerin g molecular 



cloud formation (e.g., iGir. Blitz fc Magnanilll994j ; lYamamoto et al.l 120031 . 120061 ). However 



there is still a missing link between the neutral medium and molecular clouds observed with 
the CO emission. Even if significant H_2 molecules are formed in the shell, the intermediate 
cold gas phase without CO molecules is observationally "dark". Some recent studies have 
repor ted signs of such "dark" gas phases (lOnishi et al.l l200ll ; iGrenier. Casandjian fc Terrier 
2005h . 



In our previous papers (IHosokawa &: Inutsukall2005l . l2006al . hereafter, Papers I and II), 
we have studied the expanding H II regions in the molecular clouds focusing on the detailed 
structure of the swept-up shell. We have shown that most of the swept-up gas remains in 
the shell as cold and dense molecular gas, and star-formation triggering is expected with 
various central stars and ambient number densities. After the destruction of the cloud, on 
the other hand, the stellar UV/FUV radiation spreads to the diffuse neutral medium. In 
this paper, we explore how the stellar UV/FUV radiation works for the cycle of the ISM. We 
calculate the time evolution of the shell as well as of the H II region in the ambient neutral 
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medium. While the expanding H II region photoionizes the neutral medium, the swept-up 
neutral medium may be recycled into the molecular gas. We calculate models with various 
ambient number densities (CNM and WNM), and central stars to explore the role of the 
H II regions expanding in the diffuse ISM. We show that the cold (T ~ 30 — 40 K) and 
dense (n ~ 10 3 cm -3 ) shell is sometimes formed around the H II region, and H2 molecules 
accumulate there without producing any CO molecules. This is just the intermediate gas 
phase between the neutral medium and the molecular clouds. The CO emission is not 
available to enable the search of this cold HI/H 2 gas. Instead, we try to hunt the predicted 
"dark" hydrogen with the HI self-absorption (HIS A). We find the shell-like HISA around the 
giant H II regions, W4 and W5. Our quantitative analysis shows that this HISA cloud can 
be the intermediate gas phase, just as predicted in our numerical calculations. 

The organization of this paper is as follows: In § [21 we briefly review our numerical 
method. The calculated models are also listed in this section. In § [3j we present the results 
of the numerical calculations. In § 13.11 and 13.21 we examine the expansion of the H II region 
in the ambient CNM and WNM respectively. The findings based on the observational data 
are shown in § HI § [5] and [6] are assigned to the discussions and conclusions. 



2. Numerical Modeling 



We solve the UV/FUV radiative transfer and thermal/chemical processes in a time- 
dependent hydrodynamics code. We have described our numerical code in detail in Paper II, 
and only mention some points here. The numerical sche me of the hydro dynamics is based 
on the second-order Lagrangian Godunov method (e.g., Ivan Leerl 119791 ). The frequency- 
dependent UV/FUV radiative transfer is solved with on-the-spot approximation. We include 
17 heating/cooling processes in the energy equation; for example, H-photoionization heat- 
ing, photoelectric heating, H-recombination cooling, and line cooling via Ly-a, 01 (63.0/xm, 
63.1/xm), Oil (37.29/im), CII (23.26/xm, 157.7/im), rovibrational transitions of H 2 and CO 
molecules and so on. The energy transfer by thermal conduction is also i ncluded in the energy 
equa tion to resolve the Field length of the thermally unstable media (IKoyama fc Inutsuka 
20041 ). The minimal set of the chemical reactions is implicitly solved for e, H, H + , H2, C + , 
and CO species. The si mplified reaction scheme between C + and CO is adopted, following 
Nelson fc Langerl (119971 ). Since the size distribution of grains is uncertain, we include the 
dust absorption only outside the H II region. We have confirmed that dust absorption in 
the H II region is not imp ortant in the diffuse ambient medium, o wing to the low column 
density of the H II region fjPetrosian et al.lll972l ; lArthur et al.l 120041 ) . 



We add the effect of the FUV background radiation field to our code. The gas tern- 
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perature and the ionization rate in the ambient neutral medium are initially set at the 
equilibrium v alues for the given density and FUV background field normalized in Habing 
units, Gfuv (IHabind Il968l ). If a dense shell forms around the H II region, molecular gas 
may form from the swept-up neutral gas. However, the reformed molecules in the shell can 
be photodissociated by FUV photons, both from the central star and from the background 
field. The penetration of the FUV background radiation into the s hell is treated by the an- 



alytic shielding function, measuring the optical depth from the SF (jDraine fc Bertoldilll996 



Lee et al.lll996l ). Some recent codes for the steady PDR have included more detailed calcu- 



lation of the level populations and line transfers, and shown that the resultant dissociation 
rates pos sibly deviate from those obtained by th e approximate treatment of the shielding ef- 



fect (e.g., [Shaw et al.ll2005l ; iLe Petit et al.ll2006l ). For comparison, we solved the structure of 



the irradiat ed slab that resemble s the shell appearing in our calculations, using the Meudon 
PDR code (ILe Petit et al.ll2006l ). We confirmed that these accurate calculations also show 
the physical/chemical structure very similar to that of the shell explained below. 

We adopt Gfuv = 1.0 as the standard value of the FUV background field, and include 
no attenuation in the ambient mediu m. The equilibrium n -P and n-T curves for (jfu v = 1.0 
are presented in Figure [T] (also see IWolfire et al.l Il995l ; iKoyama fc Inutsukal |2000| ). The 
cooling rate exceeds the heating rate above the equilibrium curve, and vice versa. The 
thermally unstable equilibrium states for Gfuv = 1-0 are denoted by the dotted lines. As 
the upper panel of Figure [T] shows, the two thermally stable equilibrium states can coexist 
under the same gas pressure 

of 1.5 < P < 3.2 (10- 13 dynes/cm 2 ), which are the CNM 
and WNM. The thermal balance is achieved mainly between the photoelectric heating and 
radiative cooling by [C II] 157.7/im (Ly-a and [O I] 63.1/iin) in the CNM (WNM). In Figure 
[TJ we also plot equilibrium curves with the stronger FUV fields of Gfuv = 3.3 and 10. 
The equilibrium temperature increases with the stronger FUV field, owing to the efficient 
photoelectric heating. In our numerical models presented in § El the FUV radiation field 
can be as high as Gfuv > 10 just outside the H II region, owing to the proximal massive 
star. As shown in § 13.11 below, the SF emerges and a dense shell forms around the H II 
region expanding in the ambient CNM. These equilibrium curves are useful to interpret the 
calculated physical states in the shell. We analyze the transition of physical states across 
the SF with these equilibrium curves in § 13.11 

We calculate the time evolution of the expanding H II regions in both the CNM and 
WNM. For the ambient CNM and WNM, we adopt the typical number densities of nn,o = 
10 cm -3 and n^o = 0.1 cm" 3 respectively. The equilibrium temperature and ionization 
fraction are 106 K and 5.9 x 10~ 4 in the CNM, and 7900 K and 2.9 x 10" 2 in the WNM. 
The time evolution is calculated until 10 Myr in all models. In order to focus on the role 
of stellar UV and FUV radiation, we do not include other kinetic effects produced by wind- 
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Table 1. Model Parameters 



Model 


M* (M Q ) a 


n H ,o(cm" 3 ) b 


logS uv (s- 1 ) c 


logSpuvts- 1 ^ 


Rst (pc) e 


t dyn (Myr)' 


CNM-S64 


63.8 


10 


49.35 


49.16 


18.16 


1.67 


CNM-S41 


40.9 


10 


48.78 


48.76 


11.73 


1.08 


CNM-S22 


21.9 


10 


48.10 


48.33 


6.96 


0.64 


CNM-S18 


17.5 


10 


47.02 


47.99 


3.04 


0.28 


CNM-S12 


11.7 


10 


45.38 


47.39 


0.86 


0.08 


WNM-S41 


40.9 


0.1 


48.78 


48.76 


252.7 


23.27 


WNM-S18 


17.5 


0.1 


47.02 


47.99 


65.45 


6.03 


WNM-S12 


11.7 


0.1 


45.38 


47.39 


18.59 


1.71 



a : mass of the central star, b : number density of the ambient gas, c : UV-photon number 
luminosity of the central star, d : FUV-photon number luminosity of the central star, e : 
initial Stromgren radius, f : dynamical time, R st /Cn, where Cu is the sound speed at a 
temperature of T = 10 4 K. 



driven bubbles and supernova explosions. We briefly discuss these other dynamical processes 
in §0 



3. Results of the Numerical Calculations 
3.1. Expansion in Cold Neutral Medium 

3.1.1. Fiducial Model : Expansion around 41 M & Star 

First, we consider the time evolution around a massive star of 41 M in the ambient 
CNM (model CNM-S41) as the fiducial model. Figure [2] shows the hydrodynamical evolution 
of the H II region and shell in this model. The basic time evolution is the same as that of 
the H II region in the ambient molecular material (Papers I and II). When the IF passes 
the initial Stromgren radius (~ 10 pc), the SF emerges in front of the IF. Since the gas 
temperature of the ambient CNM is only ~ 100 K, the H II region (T ~ 10 4 K) has sufficient 
excess pressure to drive the SF. The SF sweeps up the ambient CNM, and a dense shell is 
formed just around the H II region. The time evolution of the expansion is well approximated 
by, 




R lF (t) = R st M+-W- . (i) 



- 6- 



The characteristic length scale, R st is the initial Stromgren radius, 



-S'uv A 1/3 ( no \ ~ 2/3 



and timescale, tdyn is the dynamical time, 

^JWCb^M^^J , (3) 

where Suv is the UV-photon number luminosity of the central star, Cu is the sound speed 
in the H II region, and n is the ambient number density. The gas density in the shell is 
about 100 times as dense as the ambient CNM. The gas temperature in the shell is only 
~ 30 K, which is much lower than the ambient temperature. The ambient neutral medium 
just outside the H II region is slightly warmer than that in the outermost region. This is due 
to photoelectric heating by the intense FUV radiation from the central star. For example, 
the FUV flux just outside the SF is Gfuv ~ 10 at t ~ 1 Myr. Because of the geometrical 
dilution, the FUV field outside the SF gradually decreases as the HII region expands. 

Figure [3] shows the time evolution of the chemical structure of the shell. The hydrogen 
molecules gradually accumulate in the shell, owing to the rapid reformation on the surface 
of dust grains. About 45% of hydrogen atoms within the shell exist as H 2 molecules at 
t ~ 4.5 Myr, which is about the lifetime of the adopted main-sequence star. On the other 
hand, carbon monoxide molecules hardly form in the shell. The reason for this difference 
is that in this model, the column density of the shell is only N s ^ ~ 2.5 — 5 x 10 20 cm -2 , 
which corresponds to Ay ~ 0.125 — 0.25. This means that the incident FUV radiation is 
not significantly attenuated by the dust absorption in the shell. Since the dust absorption 
does not work, the reformed molecules are protected against FUV photons only by the self- 
shielding effect. In the equilibrium photodissociation region (PDR) where the self-shielding 
is dominant over the dust absorption, the column density of the hydrogen nucleon between 
the IF and H2 dissociation front (DF) is, 

iV DF ~5xl0 20 cm- 2 (^) 4/3 , (4) 



where \ = Gfuv/^ (IHollenbach &: Tielenslll999l ). For the accumulation of H2 molecules in 
the shell, N^p must be smaller than the column density of the shell, which corresponds to 
X < 2 — 4 x 10~ 2 . In our calculation, x is always only a few x 10~ 3 at the IF (see also Figure 
EJ), and the reformed U2 molecules are reliably self-shielded in the shell. In addition to the 
sufficient self- shielding in the equilibrium state, the H2 reformation timescale is sufficiently 
short in the shell. The reformation time is given by 



n 



t form ~100Myr ( — - - ) * (5) 
M0 cm -1 ' 
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(e.g., iHollenbach et al.lll97ll ). and £f orm ~ 1 Myr in the shell. Therefore, equilibrium H2 
abundance is achieved during the expansion. On the other hand, the self-shielding of CO 
molecules is inefficient owing to the small abundance of carbon and oxygen atoms. Even 
if H 2 molecules can accumulate in the shell, CO molecules cannot form without sufficient 



iV S h for shielding FUV photons by the dust absorption (jBergin et al.ll2004j . Paper II). The 
shielding by H 2 molecules is an effect secondary to the dust shielding. These features do not 
change in the H I clouds with nn,o ~ a few xlO cm -3 , but change in the molecular clouds, 
where the ambient number density is high, nn.o > 100 cm -3 . This is well understood with 
the n-dependence of the column density of the shell, iV s h and the incident FUV flux at the 
IF, F FUVii at a given t/t dyn , 

N sh oc nR st oc n 1/3 , (6) 
F FUV ,i oc Spuv/^st oc n A/3 . (7) 

As equation flBJ) shows, the column density of the shell rises with the higher ambient number 
density. Although equation (|7j) shows that the incident FUV flux also increases with in- 
creasing n, the intense FUV radiation can be shielded. Since the dust attenuation works as 
exp(— Td), the dust shielding effect quickly becomes significant once the dust opacity of the 
shell, Td (oc Agh) exceeds unity. When the H II region sufficiently expands in the molecular 
cloud, therefore, the FUV radiation can be blocked by the dust absorption, and both H 2 and 
CO molecules can accumulate in the shell (see Papers I and II). 

In our calculation, the swept-up CNM (T ~ 100 K) is converted to the cold (T ~ 30 K) 
partially atomic and molecular gas behind the SF. The formation mechanism of such a dense 
and cold shell is not understood within the framework of the isothermal SF. Furthermore, 
Figure [2] shows that the shell density and temperature gradually decrease as the H II region 
expands. We can understand all these features with the equilibrium n-P and n-T curves of 
the neutral medium. In Figure HJ "o" signs represent the physical states in the pre-shock 
region at t = 1,3, and 9 Myr in our numerical simulation. This figure shows that "o" sign 
vertically goes down with increasing time t, which means that the pre-shock temperature and 
pressure decrease almost isochorically. This is because the FUV radiation field just outside 
the shell gradually diminishes owing to the geometrical dilution, so that the photoelectric 
heating becomes inefficient. Since the heating timescale is much shorter than the dynamical 
timescale, tdyn ~ 1 Myr, the thermal equilibrium is achieved at each time step. In the same 
figure, "•" signs represent physical states in the post-shock region at each snapshot. The 
upper panel shows that the gas pressure increases by about ten times across the SF. We 
can understand this high post-shock pressure by recalling the derivation of equation flTJ. 
Equation flTJ is derived from the equation of motion of the shell, 

M sh R W = AnRUPu - ^rarn) (8) 
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(Paper II). In this equation, M s h is the mass of the swept-up shell 

Air 



M sh = yPoi?i F , (9) 



Pu is the thermal pressure in the H II region, 



3/2 



Pu = PaCfj ^g-J , (10) 
P ram is the ram pressure, 

Pram = P0-RlF> (H) 

and po is the ambient mass density. As equation ([8]) shows, the motion of the shell is 
controlled by the balance between the thermal and ram pressure from the inner and outer 
edge of the shell. Using equations ((H), ( fTUl) and ( |TTi) . the ram pressure can be written as, 

4 ( 7 [4 t 

Pram = ^p C u 1 + ~\ ~~ 1 = ^PlI- (12) 



3 \ 4 V 3t d 



yn 




The backward ram pressure is always greater than the forward thermal pressure by a factor 
of 4/3, and the expansion is decelerating. Since the sound crossing time of the shell is much 
shorter than t dyn , the pressure within the shell smoothly connects with P ram (P n ) at the outer 
(inner) edge of the shell. Therefore, the post-shock pressure is almost equal to P ram at each 
time step. The pressurized neutral medium is not initially in the thermal equilibrium state, 
but going to the equilibrium state over the cooling timescale, t coo i. Since t coo i ~ 0.1 Myr is 
much shorter than t^ym the gas element engulfed in the shell quickly reaches the thermal 
equilibrium state. Figure 0] shows that the calculated density and temperature in the post- 
shock region are equal to the equilibrium values for the given post-shock pressure (~ P ram ) 
and FUV field at each time step. 

While the gas element is going to the new equilibrium state, the total cooling rate 
significantly exceeds the heating rate. The dominant cooling process for this non-equilibrium 
gas is the line cooling of CII (157.7pm). The frequency-integrated emergent intensity of this 
line emission is about, 

~ l P ° V * = 03 X 10 " {wt/l) 3 ^/™ 2 /see/ 8t r (13) 



(IHollenbach fc McKed Il979l ). where po is the ambient mass density, and v s h is the shock 
velocity. As equation ( T131) shows, the emission from the non-equilibrium cooling layer quickly 
diminishes as the shell decelerates. Besides the non-equilibrium cooling layer just behind 
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the SF, most of the gas within the shell is in the thermal equilibrium state, where the 
photoelectric heating is mainly balanced by the CII (157.7/xm) line cooling. Therefore, the 
CII line emission also comes from the shell itself. The emergent intensity from the equilibrium 
shell is about, 



J e ~ — A C iiiV sh 

Z7T 



1.6 x 10" 



A 



en 



10 25 ergs/sec / \10 20 cm 



sli 



ergs/cm 2 /sec /str, (14) 



where Acn is the average CII cooling rate in the shell. As expected from equations (113]) and 
( I14p . the total CII emission is dominated by that from the shell in the thermal equilibrium, 
but for the very early phase of the expansion. Although this far-IR CII emission could be 
detected, this line, an ordinary feature of the PDR, is not a good signpost of the molecular 
cloud formation. Instead, we seek the predicted gas phase using the HISA feature in § HI 

In our 1-D calculation, the neutral medium is steadily swept up into the shell. However, 
some instabilities actually deform the shell and excite its internal structure. First, the dense 
gas layer is expected to fragment by the gravitational instability. The shell-fragmentation oc- 
curs first at the length scale of t he layer thickness, with a growth timescale of ~ (Gp) -1 ^ 2 (e.g.. 



Elme green fc Elmegreenlll978l ; iMiyama. Narita &: Hayashilll987l ; iNagai. Inutsuka &: Miyama 
19981 ). We adopt t > (Gp)~ 1 / 2 as the simplified criterion of the gravitationally unstable 
layer. Figure [3] shows that the shell becomes gravitationally unstable soon after t ~ 2 Myr, 
when H2 molecules begin to accumulate in the shell. The typical radius and mass of the 
fragment are about 0.1 pc and 10 M & respectively. If each fragment subsequently con- 
tracts, the molecular reformation will be promoted by the self- shielding of cores, owing to 
the increased column density. Second, the compressed layer collapsing due to the significant 
radia t ive co oling just behind the SF, suffers from thermal instability (IKoyama &: Inutsuka 
20001 . 120021 ). In Figure HJ we plot the evolution track of a fluid element across the SF 
at t ~ 3 Myr. The gas temperature rises up to T ~ 2000 K at the SF, due to com- 
pressional heating, and falls to T ~ 30 K, due to the significant radiative cooling behind 
the SF. The timescale to settle into the new equilibrium state is about the cooling time, 
^cooi ~ 0.1 Myr. The upper panel of Figure @] shows that the collapsing layer evolves al- 
most isobari cally. The linear sta b ility a nalysis of the isobarically collapsing layer has been 
explored by iKoyama fc Inutsukal (120001 ) , and the shaded region in the lower panel of Fig- 
ure H] represents the expected thermally unstable region with Gruv = 3.3. Note that the 
gas is cooling and isobarically contracting in the unperturbed state in their linear analysis. 
Whether the small perturbation grows or not is analy zed in the evolving background state, 
which differs from the original analysis by IFieldl (119651 ). Figure @] shows that the compressed 
CNM becomes thermally unstable on its way to the new equilibrium state. The linear anal- 
ysis shows that the layer f ragments into tiny cloudlet s, whose size is about tens of AUs, over 
the timescale of 10 3-4 yr (IKoyama &: Inutsukal l2000l ). Mult i- dimensional numerical simula- 
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tions (e.g., iKoyama Inutsuka] 120021 ; lAudit &: Hennebelld 120051 ) show that each fragment 



has some translational motion of a few km/s. This motion is driven by the pressure of the 
surrounding warmer, less dense region. Thus, the motion is inevitably subsonic for the sur- 
round ing medium, and does not dissipate over the sound- crossing time. IKoyama &: Inutsuka 
( 120021 ) have suggested that this translational cloudlets motion is the origin of the observed 
broad line width, or "turbulence" in the cold clouds. Our current 1-D calculations cannot 
explain the the actual evolution of the fragmentation of the shell, which will be affected by 
both the thermal and gravitational instabilities. Future multi-dimensional simulations of 
radiative-hydrodynamics will reveal the detailed nature of the fragmented shell. 



3.1.2. Dependence on Central Stars 

Second, we study how time evolution changes with different central stars in the same 
ambient CNM. Our numerical results show a clear feature; with the higher-mass central star, 
the shell density is higher and the reformation of H 2 molecules becomes more efficient in the 
shell. 

For example, Figure [5] presents gas-dynamical evolution around a massive star of M* = 
11.7 M & (model CNM-S12). The SF emerges and sweeps up the ambient CNM, which is the 
same as model CNM-S41. The shell density is 10 times as high as the ambient CNM density 
at t ~ 1 Myr, but quickly declines afterward. The shell density is always lower than that in 
model CNM-S41 at the same time, t (c.f. Figure [7j). Despite the different densities of the 
shells, the temperature profiles in Figures [2] and [5] are similar. This is because the FUV field 
outside the shell is similar among these models (Figur^TJ). The FUV flux just outside the SF 
is Gfuv ~ 10 at t ~ 1 Myr, and gradually decreases owing to the geometrical dilution. The 
formation of molecules in the shell hardly occurs in model CNM-S12. No CO molecules are 
formed, and only less than 0.1% of H atoms is converted to H2 molecules in the shell over 
10 Myr. Since the column density of the shell is A" s h < 2.5 x 10 20 cm -2 , the dust absorption 
does not protect reformed molecules against FUV photons. As shown in equation (JH), the 
self- shielding becomes less efficient with the larger x — Gfuv/^ in the shell. In model CNM- 
S12, x is about ten times as large as that in model CNM-S41 at the same time, t (also see 
Figure [7]). Therefore, A^df given by equation (0J exceeds the column density of the shell, and 
even the self-shielding of H 2 molecules does not work. Furthermore, if the self-shielding were 
efficient, H 2 molecules would not form in the low-density shell. Equation §5§ shows that 
the reformation timescale becomes longer than 10 Myr, when the shell density is smaller 
than 100 cm -3 . The reformation timescale is so long that no equilibrium H2 abundance is 
achieved. We have calculated other models with different exciting stars, and found that the 
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formation of H2 molecules is more efficient with higher-mass stars. Figure [6] shows the time 
evolution of the chemical structure of the shell around the central star of 17.5 M (model 
CNM-S18). In this model, H 2 molecules are slightly formed in the shell. About 6 % of 
swept-up hydrogen atoms are converted to H 2 molecules by the time of 10 Myr. 

We investigate why the efficiency of H2 formation depends on the mass, or UV/FUV 
luminosity of the central star. The key physical quantities are the ratio, x = Gp\jy/n and 
number density, n, in the shell. The former reflects the efficiency of the self-shielding effect, 
which determines the H 2 abundance in the equilibrium state, and the latter determines the 
reformation timescale, over which the H2 abundance becomes the equilibrium value. Figure 
[7J shows the time evolution of the number density, FUV radiation field, and their ratio, x i n 
the shell in models CNM-S41, S18, and S12, respectively. At the same time, t, the number 
density is larger with higher-mass stars, and the FUV fields are similar among models. Thus, 
their ratio, x is smaller with a higher-mass star. The self-shielding effect is more efficient 
with a higher-mass central star, as equation (j3J) shows. The higher shell density ensures that 
the equilibrium state is achieved over a shorter timescale, and that H2 molecules quickly 
accumulate in the shell with a higher-mass star. Consequently, the variation in the chemical 
structure of the shell originates in similar and different time evolutions of Gfuv an d n among 
models shown in Figured 

We make use of some analytic relations to interpret these characteristic features. First, 
the FUV field at time t, Gfuv(^) scales with the UV/FUV luminosity of the star, 

Gpuv(()K *- K ^(^)-( 1+ z^y 8/7 . (15) 

This shows that Gfuv(^) is normalized with the first term, Sp^ v ( S'fuv /S\jv) 2 ^ 3 , which is 
smaller with the more massive star. On the other hand, the last evolutionary term has an 
inverse dependence, because the dynamical time, td yn is longer with the higher-mass star (see 
equation Q). Owing to the cancellation of these opposite dependences, the time evolution 
of the FUV field is similar among models with different stars, as shown in Figure [7J 

Next, let us consider why the shell density is larger with a higher-mass star. In § I3.1.1[ 
we have shown that the post-shock density is determined by the FUV radiation field and 
gas pressure in the shell at each time step. Because of the similarity in the time evolution of 
Gfuv among models, only gas pressure in the shell dominates the post-shock physical state. 
The time evolution of the ram pressure is given by equation f)12p . where the dynamical time, 
tdyn is longer with the higher-mass star. Therefore, thermal pressure in the shell is higher 
with the higher-mass star at the same time, t, which causes a variation in shell density 
among models. Figure [8] presents the time evolution of thermal pressure in both pre-shock 
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and post-shock regions in models CNM-S41, S18, and S12. Since the FUV radiation field is 
similar, the time evolution of the pre-shock pressure is almost the same among models. The 
post-shock pressure is higher with the higher- mass star, as predicted by equation f[T2l . The 
significant enhancement of the pressure behind the SF accounts for the high density of the 
shell. Figure [9] presents several post-shock and pre-shock physical states in model CNM-S12. 
In this model, the post-shock pressure is just a little higher than the pre-shock pressure, and 
the shell density is significantly lower than that in model CNM-S41 (c.f. Fig. H]). 

Finally, we discuss the multi-dimensional structure of the shell expected to be produced 
by various instabilities. First, gravitational instability grows earlier in denser layers. Since 
the shell density is higher with a higher-mass star, shell-fragmentation will occur earlier 
in these higher-mass models. While the shell becomes unstable at t ~ 2 Myr in model 
CNM-S41 (Figure [3]), it is not until t ~ 5 Myr that the unstable region appears in model 
CNM-S18 (Figure [6]). No gravitationally unstable region appears in CNM-S12. Second, the 
thermal stability of the isobaric cooling layer depends on the compression level across the 
SF. Figure [9] shows the evolution track of the layer at t ~ 3 Myr in model CNM-S12. The 
pressure rises by a factor of two or three across the SF, and the evolution track only grazes 
the thermally unstable region on the n-T plane. In summary, our calculations show that 
H 2 formation is more efficient in a model with a higher-mass star, and that the effects of 
gravitational/thermal instability are more significant in such higher-mass models. 



3.2. Expansion in Warm Neutral Medium 

Once the stellar UV radiation escapes from the molecular clouds and H I clouds (CNM), 
the H II region expands to the WNM. Figure [TD] shows the hydrodynamical evolution of the 
H II region around a 41 M & star in the ambient WNM (model WNM-S41). Because of 
the low ambient number density, the H II region expands into a very large region of more 
than 100 pc. As Figure [10] shows, however, the expanding H II region hardly affects the 
hydrodynamics. The gas density slightly increases by a factor of two in front of the IF. 
This is partly because of the large initial Stromgren radius. In model WNM-S41, the initial 
Stromgren radius is, R st ~ 250 pc, and the H II region remains in its early formation phase 
at t < 2 Myr. Furthermore, after the H II region reaches the initial Stromgren radius, no 
dense shell is formed around the H II region. The ambient gas temperature is originally as 
high as T ~ 8000 K, and the pressure excess of the H II region over the ambient WNM is 
small. No SF with a high Mach number ever emerges in front of the IF. The ambient WNM 
is only converted to the ionized medium with a similar temperature and density, which is 
the WIM. These basic features of the time evolution do not change with different central 
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stars. Although iHennebelle fc Peraultl (119991 ) and iKoyama Inutsukal (120001 ) have shown 
that the highly pressurized WNM can be converted into the CNM, this is not the case with 
the expanding HII regions. The bubble driven by stellar-wind or supernova explosion can 
provide the much higher pressure even in the WNM, and triggers the phase transition from 
WNM to CNM (also see § El- 



4. Hunt of the Dark Molecular Hydrogen Clouds in the Observational Data 
4.1. Shell- like HI Self-absorption around Giant H II Regions 



When the H II region expands in the CNM, the SF emerges in front of the IF and 
the ambient CNM is compressed in the shell. We have shown that H 2 molecules can form 
if the central star is more massive than 18 M Q , but CO molecules seldom form in the 
shell. The calculated temperature and density in the shell are T ~ a few x 10 K and 
n ~ 100 — 1000 cm -3 . This is just the gas phase intermediate between the ambient CNM 
and molecular clouds observed with the CO emission. However, it is generally difficult to 
detect such an intermediate gas phase. The gas temperature is too low to thermally excite the 
rovibrational transition of H 2 molecules. The FUV radiation field in the shell is Gfuv ~ 10 
at most, and it is hard to detect the UV fluorescent H2 emission from the distant giant 
H II regions. Instead of hunting for signs of the cold H2 molecules, we probe the predicted 
"dark" hydrogen with the H I self- absorption (HIS A). The HIS A is the H I 21-cm line- 
absorption against the background emission from the warmer neutral medium. Statistical 
studies of the CNM usually analyze the distribution of the 21-cm absorption features against 
the st rong radio continuum sources, such as comp act H II regions and supernova remnants 
(e.g., bickev et all I2OO3I : iHeiles fc Trolandl l2003al lbh. On the other hand, the brightness 
temperature of the background H I 21-cm emission is only Tb ~ 100 K, and HISA features 
automatically trace the colder tail of the CNM. 

We have sought the HISA features associated with the giant H II regions using recen t 
observational data by the Canadian Galactic Plane Survey (CPGS; see iTaylor et al.l 120031 ). 
The CGPS data provides us with the high-resolution (~ 1') maps of HII, HI, CO, and dust in 
the Galactic plane. We use the 21-cm continuum image and velocity cube of the 21-cm emis- 
sion by the Dominion Radio Astrophysical Observatory (DRAO) survey (IHiggs fc Tapping 
2000h . 12 CO(1-0) data cubes by the Fiv e College Radio Astronomy Observatory (FCRAO) 
outer g alactic survey (Heyer et al.lll9981 ). and reprocessed IRAS images at 60 and 100 fim 
bands (ICao et al.l 119971 ). We have examined the star-formation complex W3-W4-W5 along 
the Perseus arm, which includes the gian t H II regions, W3, W4 and W5. Adopting a 
distance of 1.95 kpc to the Perseus arm (IXu et al.l 120061 ). the angle of 1° corresponds to 
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32.2 pc. Figure [TT] shows the multi-wavelength view of the W5 region. The upper left panel 
is the 21 cm continuum map, which shows distribution of ionized hydrogen. The upper-right 
panel presents a 60 /im image of the W5 region, and we can see the clear shell-like struc- 
ture around the ionize d gas. Other IRAS-baad images also show similar shell-like structures 



( iKarr fc Martini 120031 ) . The lower left panel of Figure [TT] shows a channel map of 21-cm 
emission at ^lsr = —39.8 km/s. We have also found a narrow shell-like feature in this 
map, which shows strong spatial correlation with the dust shell, as indicated by the white 
rectangle in the lower right panel. The brightness temperature along the shell-like feature is 
Tb ~ 70 K, which is lower than that on both sides of the shell, Tb ~ 100 K. Since it is very 
unlikely that the 21-cm emission from the diffuse WNM would fluctuate over such a small 
scale, and that such fluctuation would overlap with the dust shell by accident, we consider 
this shell-like feature as being HISA by the colder neutral medium. Furthermore, these dust 
shell and shell-like HISA feature should originate in the same structure just around the W5 
H II region. We also show the contours of the 12 CO(1-0) channel map at the same velocity 
in the upper right panel, but the CO emission does not trace the shell-like structures with 
the current sensitivity of the FCRAO survey. 

We have also found some similar features in the giant H II region, W4. Figure [T21 shows 
a multi-wavelength view of the W3-W4 region. The upper left panel shows the distribution 
of the ionized hydrogen. The W4 H II region has an ionized shell-like structure. The W3 
H II region is denser and more compact than W4. The upper right panel presents the 60 /im 
image superposed on the contours of the 12 CO(1-0) channel map at t> L sR = —47.3 km/s. 
The distributions of both the dust and CO emission are two-fold. In I > 134.3°, the 60 /im 
emission shows a shell-like structure just outside the ionized shell. The CO molecules are also 
along the dust shell, but discretely distributed as small clouds. In I < 134.3°, on the other 
hand, the broad and strong 60 /im emission comes from the warmer dust around the W3 H II 
region. The CO emission also broadly extends in I < 134.3°. This is the eastern edge of the 
W3 giant molecular cloud (GMC), which actually extends over / < 133°. The lower-left panel 
of Figure [T2l shows the 21-cm channel map at t>LSR = —47.3 km/s. The HISA feature has also 



been f ound, and overlaps with the dust shell (lower- right panel). ICarpenter. Heyer fc Snell 



( 120001 ) have shown that the small molecular clouds are widely distributed in the W3-W4-W5 
cloud complex, and have speculated that these are remnants of pre-existing molecular clouds 
dispersed by the negative feedback from massive stars. We suggest that a number of the 
small clouds form from the compressed ambient neutral medium. 
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4.2. Quantitative Analysis of the Data 

In order to constrain the physical state of the observed HISA cloud, we quantitatively 
analyze the data. We define "on"- and "off" -regions designated by 6' x 6' squares in the 
lower-left panel of Figure [ITJ The on-region is just on the HISA feature around the W5 
H II region, and the off-region is slightly outside the HISA feature. Figure d3] presents the 
averaged brightness temperatures in the on-region (T on ) and off-region (T q) at different 
velocities. The emission profile in the on-region clearly shows a depression compared with 
that in the off-region. We regard this absorption feature as the HISA profile, and obtain 
the line profile calculating T on — T oS . The line depth and width are about AT ~ 30 K, 
and Av ~ 5 km/s respectively. The velocity at the line center is vlsr ~ — 40 km/s, which 
is the velocity of the channel maps presented in Figure [TlJ The HISA profile is available 
to limit This a and thtsa consid ering the radiative transfer under a given geometry of the 
system (e.g.. iGibson et al.ll2000l ). We suppose that the HISA cloud, whose spin temperature 
and optical depth are T H isa and t H isa, lies in the sea of the WNM. The diffuse WNM is 
distributed in both the foreground and background of the HISA cloud. The spin temperature 
and the optical depth of the foreground (background) WNM are represented as Tf g (T bg ) and 
r fg ( r bg)- Note that the spin temperature of the WNM differs from the kinetic temperature 
by a factor of two or three, be cause the t hermal equilibrium is not achieved owing to its 



low density, n ~ 0.1 cm 3 (e.g.. lLisztll200ll ). The continuum radiation includes the free-free 



emission from the H II region, cosmic microwave background and so on. In our case, the W5 
H II region is so diffuse that the brightness temperature of the continuum radiation, T c , is 
much smaller than that of the line emission. Presuming that all the gas components except 
the HISA cloud are optically thin, the radiative transfer equation along the line of sight in 
the on- and off-region are, 



T on — i~f K T {R + T H isa(1 — e 



THISA 



+ r bE T be e-™ SA + T c e~~ miSA - T c 



(16) 



T Q g — T fg T fg + T bg T b g, (17) 

where T on (T ff) is the brightness temperature observed in the on- (off-) region. Subtracting 
equation (JTTj) from equation ffTBI) . we obtain, 



T n T, 



off 



;i-e- THISA )(T H isA-pT off -T c 



(18) 



Gibson et allboool : iKavars et al1l2003l : Eerto^SoJ, where p is the ratio of the back- 



log., 

ground WNM emission to the total emission, p = T Dg T hg /T o{i . We derive the relation between 
^hisa and thisa fc> r a given p from equation f[T8l . because T on , T oS and T c are all observable 
quantities. We adopt the numerical values of T on = 67.9 K, T off = 97.3 K at the line center, 
and T c = 9.1 K. In Figure [HJ the solid lines represent the Thisa-Thisa relations with these 
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values for a given p. As this figure shows, we can set the maximum spin temperature of the 
HIS A cloud as T H isa < 77 K, where the equality is satisfied when p — 1. Inversely, if the fore- 
ground WNM emission accounts for more than 25% (50%) of the total WNM emission, which 
corresponds to p < 0.75 (0.5), the spin temperature should be lower than 52.7 K (28.3 K). 
The estimated temperature is comparable to the temperature p redicted by our calculations 
and much lower than that of the typical PDR, T ~ 100 — 1000 K (jHollenbach fc Tielendll999 
Papers I and II). 



We can set further constraint on the This a- This a plane using the equation 



Thisa = 5.2 x 10 
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HI 



Aw Thisa 



(19) 



( IDickey fc Lockmanlll990l ). where Ahi is the column density of atomic hydrogen in the HIS A 
cloud. The line width, Av is about 5 km/s from the line profile presented in Figure [TBI 
Therefore, another T H isa-thisa relation is obtained by equation f!T9l evaluating N m . Taking 
into account that the HISA cloud can be partially molecular, the atomic column density 
is written as Ahi = /Ah, where Ah is the total column density of all the hydrogen nuclei 
and / is the atomic fraction in the HISA cloud. We evaluate the total column density in 
two different ways, instead of directly measuring th e atomic column density. First, we use 
the extinction map given by 



the Digitized Sky Survey I (ILasker 



Dobas 



liet 



al.l (120051 ) . They have analyzed the database of 



1994! ) using the star count technique, and completed 



an extinction map of the entire the Galactic plane. We extract the Ay map at 6' angular 
resolution around the W5 region. In order to estimate the extinction only through the 
HISA cloud, we calculate the average extinctions in the on- and off-regions, and derive 
their difference. The excess extinction in the on-region is Ay = 0.81, which corresponds to 
A H ~ 1.6 x 10 21 cm -2 . With this value of A H , we plot the T H isa-thisa relation given by 
equation (Tl~9"]) for / = 1.0, 0.1 and 0.01 in Figure [TH In § 13.1.11 and 13.1.21 our numerical 
models have predicted that the gas phase of / ~ 0.5 and T ~ a few xlO K should form 
around the H II region. Figure shows that such a gas phase is still possible under our 
constraints. Second, let us evaluate the total column densit y in another way, using the 60- 
and 100 — /im dust emissions (e.g, lArce fc Goodman! Il999i ). The intensity of the thermal 



and optically-thin dust emission is given by, 



rxB u (T d ) 



(20) 



where t\ is the dust opacity, B v is the Planck function, and T<j is the dust thermal temper- 
ature. The dust opacity is written as, 



r A = e A A H , 



(21) 
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where e\ is the dust emissivity per hydrogen atom, which is usually expr essed as e\ = en A 13 . 



Although (3 and e actually vary with the different environments (e.g., iDupac et all 120031 ) 
we adopt the typical emissivity law in the local ISM, 



e A = 1.0 x 10~ 25 ( ^ cm 2 (22) 



250 fim 



(IBoulanger et al.lll996l ). Using equations (|20l) . ( J2T1) and ( J22l) . we can derive Td and Nn from 
the intensities at 60 and 100 /im, Ieo^ m and Jioo^m- We analyze the IRAS 60 — /im and 
100 — /im images, and extract the dust emission only from the dust-shell by calculating the 
difference in intensity between on- and off-regions. The excess intensities in the on-region are 
= 89.1 MJy/str and hoo^m = 180.3 MJy/str. With these values, we obtain Td = 29.5 K 
and A^h = 9.5 x 10 20 cm -2 . The derived total column density is smaller than that from the 
extinction map by a factor of 1.7 only. In Figure dU we also plot the Thisa-thisa relation 
given by equation (fi~9l) with this value of A% for / = 1.0. Note that the 60-/im emission is 
actually contaminated with the emission from very small grains which are not in thermal 
equilibrium. Therefore, the above values of Td and iVn should be considered as the upper 
and lower limits. The future high-resolution observation at A > 100 /im will provide tight 
constraints. 



5. Discussions 

5.1. Strength of the FUV Background Field 

In our numerical modeling, we have adopted an FUV background field of Gfuv = 1-0 as 
the fiducial value. However, the F UV background field is generally variab le with time, and 



sensitive to proximal massive stars (jParravano. Hollenbach fc McKed 120031 ) . Quantitatively, 



a stronger FUV field reduces the molecular abundance in the shell. This is because the 
self- shielding effect becomes inefficient with a stronger FUV background field, namely the 
ratio, x = G-pw/n becomes larger with larger Gfuv- As expected with the equilibrium n-P 
curves (also see §[2]), the shell density decreases as Gfuv increases under the same pressure, 
and x increases. We have examined some models with different FUV background fields. 
With G FUV = 3.3 (10) in model CNM-S41, for example, only about 25 % (2 %) of hydrogen 
atoms exist as U2 molecules in the shell at t ~ 4.5 Myr. Remember that the H2 ratio is 
45 % in the same model with Gruv = 1-0 (see § I3.1.ip . The gas temperature in the shell is 
about 40 K (60 K) with Gfuv — 3.3 (10), which is still much lower than the typical CNM 
temperature. This situation may be more suitable for HISA clouds around the W4 and 
W5 H II regions, because the FUV field should be higher than the local value in a massive 
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star-forming region. Inversely, the H2 abundance ratio is slightly enhanced with a weaker 
background field. However, the FUV field in the shell does not fall below Gfuv ~ 1 owing 
to the FUV radiation from the central star. CO molecules seldom accumulate in the shell, 
and the "dark" character of the shell does not change with a weaker background field. 



5.2. Reformation Rate of Hydrogen Molecules 

In previous sections, we have seen that the expanding H II region sometimes triggers the 
formation of H2 clouds without CO from the ambient neutral medium. These H2 molecules 
are formed on dust g r ains i n the dense shell, and we included the reaction rate given by 



Hollenbach fc McKed (119791 ) in our numerical calculations. However, there is still a large 



uncertainty on the formation rate of H2 molecules. It is not sur prising that the form ation 



rate depends on the local physical state of the ISM. For example, lHabart et al.l (12004 ) have 
suggested that the formation rate should be about five times enhanced in some nearby PDRs. 
Taking the uncertainty into account, we have calculated model CNM-S41 with different 
formation rates of 5 x Rf and 1/5 x Rf, where Rf is the adopted normal reaction rate. 
Whereas about 45 % of H atoms are included in H2 molecules with Rf, H2 fraction becomes 
83 % (3 %) with 5 x Rf (1/5 x Rf) at t ~ 4.5 Myr. Since the self-shielding is the dominant 
process to shield the FUV radiation, the H 2 abundance within the shell is somewhat sensitive 
to the formation rate. In spite of different H 2 abundance in the shell, however, only less than 
0.1 % of C atoms are converted to CO molecules in all models. As mentioned in § 13. 1\ 
the dust absorption of the FUV radiation with the sufficiently large N s h is needed for the 
accumulation of CO molecules. 



5.3. Expansion of Bubbles and Superbubbles 



Although we have only focused on the expansion of the H II region in this paper, the 
strong stellar wind generally affects the dynamics around a main-sequence massive star. 
When the wind-driven bubble expands in a preexistent H II region, the forw ard SF sweeps 
up a n ionized medium and the ionized shell forms around the hot bubble (IWeaver et al. 
19771 ). The ionized shell observed in the W4 H II region (Figure IT21) may correspond to 
this feature. In the diffuse ambient medium of n < 10 cm -3 , however, the effect of the 
wind-driven bubble is not so significant as for expansion around a single star. We can 
easily show that the bubble pressure quickly decreases to the initial H II pressure, an d the 
bubble is confined to the H II region (IMcKee et al.lll984l ; iHosokawa &: Inutsukall2006bl ). The 
dynamical effect of the bubble is highly enhanced when the giant bubble expands around 
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an O B association (superbubble; e.g.. iTomisaka. Habe fc Ikeuchil Il98ll ; iMcCray fc Kafatos 
19871 ). A huge energy input is maintained over several 10 Myr by the stellar winds and 
continuous supernova explosions. The expansion is mainly driven in the later supernova 
phase, and the energy input rate is L s b ~ 6 x 10 35 iV* erg/s, where N* is the number of OB 
stars. The time evolution of the superbubble size and pressure is given by, 
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(23) 



(24) 



( iWeaver et al.lll977l : IMcCray fc Kafatoslll987l ). In § 13.21 we have shown that the H II region 
causes little dynamical effect in the ambient WNM. However, equation f[2~4"l) means that the 
superbubble pressure can be about ten times as high as the WNM pressure with Gfuv ~ 1 
(also see Figured]). The SF emerges owing to the high bubble pressure, and a dense shell 
forms around the superbubble. We can expect a shell density with an equilibrium n-P 
curve, as carried out in § 13.11 The swept-up WNM should tend forward the equilibrium 
state determined by the enhanced pressure and the FUV field at each time step. Adopting a 
post-shock pressure of 10~ 12 dynes/cm 2 and Gfuv ~ 1 in the shell, we expect that the shell 
density can be as high as several xlO 2 cm -3 , though the Mach number of the SF is quite 
low. The shell temperature should be several x 10 K. The molecular abundance in the shell 
depends on the column density of the shell, which is estimated as, 

Vs , „ s 4/5 ( t \ 3/5 
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(25) 



The expected column density is much lower than that of our numerical models presented in 
§ 13.11 The situation should be very difficult for molecular formation in the shell. Not only 
CO, but H2 molecules suffer from photodissociation by the FUV background field. Even 
with the higher ambient density of n ~ 1 cm -3 , the column density is not high enough for 
the accumulation of CO molecules. The above estimation agrees with the fact that galactic 
supershells are usually detecte d with the H I and/or dus t emission, but seldom associated 
with molecular emission (e.g., lHartmann fc Barton! 119971 ) . In some regi ons, however, the 



discrete molecular clou ds are clearly distributed along the supershell (e.g., iFukui et al.lll999 



Yamaguchi et all 120011 ). If these molecular clouds are formed from the swept-up neutral 



medium, an extra shielding effect should be necessary to prevent the photodissociation of 
CO molecules. Possible processes are shell-fragmentation and subsequent contraction of each 
fragment. If each fragment contracts to a dense core, the column density increases and CO 
molecules will be protected against FUV photons. Both dynamical and chemical evolution 
of the shell and fragments should be studied in detail in a separate paper. 
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6. Conclusions 

In this paper, we have studied the role of the expanding H II region in an ambient 
neutral medium. The time evolution of the H II region and swept-up shell has been studied 
by solving the UV/FUV radiative transfer and the thermal/chemical processes in a time- 
dependent hydrodynamics code. We have analyzed the chemical structure of the shell, and 
examined how efficiently the molecular gas can reform in the shell. First, we have studied 
the expanding H II region in the ambient CNM. We have analyzed its expansion around a 
massive star of M* = 41 M as a fiducial model. 

1. The SF emerges in front of the IF and sweeps up the ambient CNM. H 2 molecules 
reform in the shell, but CO molecules seldom form. This is due to the different effi- 
ciencies of the self-shielding effects. Even if H 2 molecules accumulate in the shell, CO 
molecules cannot form without the sufficient column density of the shell. 

2. The swept-up shell becomes cold (T ~ 30 K) and dense (n ~ 10 3 cm -3 ). This physical 
state in the shell is the equilibrium of the pressurized neutral medium with the FUV 
radiation field in the shell at each time step. The gas pressure rises to about the ram 
pressure in front of the SF, just after the gas is engulfed in the shell. 

3. The shell will fragment into small clouds due to gravitational instability. The typical 
mass- and length-scale of each fragment are ~ 10 M and ~ 0.1 pc. If each fragment 
contracts into a dense core and increases the column density, CO molecules may form. 
The cooling layer just behind the SF also suffers from the thermal instability. This 
layer will fragment into tiny (~ tens of AUs) cloudlets with the translational motion 
of a few km/s. 

The physical and chemical structure of the shell clearly depend on the mass, or UV/FUV 
luminosity of the central star. The dependences are summarized as follows: 

4. The formation of H 2 molecules is more efficient for a higher-mass central star. When 
the stellar mass is less than 18 M , neither CO nor H 2 molecules form in the shell. 

5. In the model with the higher-mass star, the shell density and pressure are higher. 
The denser shell becomes gravitationally unstable earlier, and the thermal instability 
also works efficiently in the highly pressurized shell. Therefore, the effect of these 
instabilities is more significant with the higher-mass star. 

Our numerical models have shown that cold and dense HI/H 2 gas is forms in the shell without 
CO molecules, and this formation is just the intermediate phase between the neutral medium 
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and molecular clouds. In order to probe this "dark" hydrogen, we have searched for HISA 
features associated with giant H II regions using the GCPS data. Our findings are as follows: 

6. We have found shell-like HISA features around the ionized gas in both the W5 and W4 
H II regions. The shell-like HISA shows good spatial correlation with the dust shell, 
but not with the 12 CO(1-0) emission. 

7. According to our quantitative analysis of the data, the observed HISA cloud can be 
in the cold (T ~ a few xlO K) gas phase that includes a significant amount of H 2 
molecules, just as in the numerical simulations. 

Finally, we have also calculated the time evolution of the H II region in the WNM, but 
the hydrodynamical evolution is much simpler, owing to the small pressure excess of the 
photoionized gas. When the H II region reaches the WNM, therefore, the UV radiation 
simply transforms WNM into the WIM. 
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Fig. 1. — Equilibrium gas temperature and pressure of the neutral medium for a given 
gas density and FUV radiation field. In each panel, the broken, dotted, and solid lines 
respectively represent the equilibrium curves with different FUV radiation fields, Gfuv — 10, 
3.3, and 1. The dotted curve in the upper panel denotes the thermally-unstable equilibrium 
states for Gfuv = 1- 
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Fig. 2. — Snapshots of the gas-dynamical evolution of model CNM-S41. In each panel, five 
snapshots represent the profiles at t — 1, 3, 5, 7 and 9 Myr respectively. The broken line in 
panel (3) represents the initial temperature of the ambient neutral medium. 
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Fig. 3. — Time evolution of the column density of each region in model CNM-S41. The 
vertical axis represents the column density of the hydrogen nucleon from the central star, 
and horizontal axis represents the elapsed time of the calculation. Two thick solid lines 
correspond to the positions of the SF and IF, and the region enclosed between these lines 
corresponds to the shell. Although the radius of the IF monotonically increases, the column 
density of the H II region decreases owing to the decrease of the density. Thin contour curves 
represent the chemical abundance in the shell; positions where X Ha = 2nH 2 /^H nuc = 0.3 and 
0.7 (broken), and X C o = «co/nc mc — 2-0 x 1CT 4 (dot-solid). The shaded region denotes the 
gravitationally unstable region, where t > (Gp)~ - 5 . 
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Fig. 4. — The radiative equilibrium states shown on the n-P (upper panel) and n-T (lower 
panel) planes. Solid, dot-solid, and broken lines represent the radiative equilibrium curves 
for different FUV fields, Gfuv = 1) 3.3, and 10 respectively. In each panel, "o"("»") signs 
indicate the calculated states in the pre-shock (post-shock) region at t = 1, 3 and 9 Myr in 
model CNM-S41. In the upper panel, the dotted line indicates P ram at t = 3 Myr calculated 
by equation (fl2"|) . The evolutionary track across the SF at t ~ 3 Myr is shown by thin 
solid line in each panel. In the lower panel, the shaded region denotes where the isobarically 
contracting and cooling layer is expected to be thermally unstable by the linear stability 
analysis (see text). 
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Fig. 5. — Same as Figure [5] but for model CNM-S12. 
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Fig. 6. — Same as Figure [3] but for model CNM-S18. The contours with thin broken curves 
represent the positions where Xh 2 = 0.01, 0.05 and 0.1 respectively. 
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Fig. 7. — Time evolution of the density of the shell, FUV radiation field at the SF, and 
their ratio, x- The solid, broken, and dot-solid curves in each panel represent the evolution 
of models CNM-S41, S18, and S12 respectively. 
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Fig. 8. — Time evolution of the gas pressure in the pre-shock and post-shock regions. The 
solid, broken, and dot-solid lines in each panel represent the evolution of models CNM-S41, 
S18, and S12 respectively. 
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Fig. 9. — Same as Figure H] but for model CNM-S12. 
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Fig. 10. — Snapshots of the gas-dynamical evolution in model WNM-S41. In each panel, 
three snapshots represent the profiles at t — 1, 3, 5, and 9 Myr respectively. The broken line 
and arrows in panel (3) represent the initial equilibrium temperature of the ambient warm 
neutral medium, and the position of the IF at each time step. 



-36- 




Fig. 11. — Multi- wavelength view of the W5 region using CGPS data. Upper left : image of 
the 21-cm continuum emission (H II gas). The brightness temperature is linearly scaled from 
5 K (bright) to 12.5 K (dark). Upper right : image of the 60-/iin dust emission and contours 
of 12 CO(1-0) line emission at t>LSR = —39.8 km/s. The intensity of the 60-/xm dust emission 
is linearly scaled from 10 MJy/str (bright) to 200 MJy/str (dark). The contour levels of CO 
emission are Tb = 1, 2, 5, and 10 K. Lower left : channel map of the H I 21-cm emission at 
v lsr — —39.8 km/s. The brightness temperature is linearly scaled from 45 K (bright) to 
110 K (dark). The blue (white) 6' x 6' square denotes the on- region (off-region) used for 
the analysis. Lower right : channel map of the H I 21-cm emission at wlsr — —39.8 km/s 
and contours of 60-/im dust emission. The contour levels are Ieo^ m = 75, 92.2, 143.8, 230 and 
350 MJy/str. The region enclosed by the white line denotes where the HISA feature shows 
the good spatial correlation with the dust shell. 
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Fig. 12. — Multi-wavelength view of the W3-W4 region using CGPS data. The distribution 
of panels and the scaling of images and contours are almost the same as those in Figure [TTJ 
Some differences are as follows: Upper left : The brightness temperature is linearly scaled 
from 4.9 K (bright) to 20.0 K (dark). Lower right : The contour levels of the 60-/im intensity 
of the dust emission are Ieo^m = 65,82.8, 136.3,225.3, and 350 MJy/str. The channel maps 
of 12 CO(1-0) and H I 21 cm emission are at ^lsr = —47.3 km/s. The region enclosed by the 
white line denotes where the HISA feature shows the good spatial correlation with the dust 
shell. 
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Fig. 13. — Upper panel : variation in the brightness temperature with the velocity, t>LSR m 
the on-region (broken line) and off-region (dotted line) in the W5 H II region. The positions 
of the on- and off-region are presented in the lower left panel in Figure [TTJ Lower panel : 
the absorption line profile obtained by T on — T off . 
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Fig. 14. — Constraints on the spin temperature and optical depth of the HISA cloud in the 
W5 region. The solid curves present the relation between Thisa and trisa given by equation 
( |T8l) for p (= TbgTbg/T ff ) = 1.0, 0.75, 0.5 and 0.25. The dot-solid curves represent the relation 
given by equation (PJ for / (= N m /N u ) = 1.0, 0.1, 0.01 with JV H = 1.6 x 10 21 cm" 2 , which 
is estimated from the extinction map. We also plot the same relation for / = 1.0 with 
iVn = 9.5 x 10 20 cm -2 , which is estimated from IRAS 100 /im and 60 /zm images (dotted 
curve). The shaded region denotes the allowed range of Thisa and thisa under our constraints. 



